A novel approach to simulate gene-environment interactions in complex diseases

Background: Complex diseases are multifactorial traits caused by both genetic and environmental factors. They represent the major part of human diseases and include those with largest prevalence and mortality (cancer, heart disease, obesity, etc.). Despite a large amount of information that has been collected about both genetic and environmental risk factors, there are few examples of studies on their interactions in epidemiological literature. One reason can be the incomplete knowledge of the power of statistical methods designed to search for risk factors and their interactions in these data sets. An improvement in this direction would lead to a better understanding and description of gene-environment interactions. To this aim, a possible strategy is to challenge the different statistical methods against data sets where the underlying phenomenon is completely known and fully controllable, for example simulated ones. Results: We present a mathematical approach that models gene-environment interactions. By this method it is possible to generate simulated populations having gene-environment interactions of any form, involving any number of genetic and environmental factors and also allowing non-linear interactions as epistasis. In particular, we implemented a simple version of this model in a Gene-Environment iNteraction Simulator (GENS), a tool designed to simulate case-control data sets where a one gene-one environment interaction influences the disease risk. The main aim has been to allow the input of population characteristics by using standard epidemiological measures and to implement constraints to make the simulator behaviour biologically meaningful. Conclusions: By the multi-logistic model implemented in GENS it is possible to simulate case-control samples of complex disease where gene-environment interactions influence the disease risk. The user has full control of the main characteristics of the simulated population and a Monte Carlo process allows random variability. A knowledge-based approach reduces the complexity of the mathematical model by using reasonable biological constraints and makes the simulation more understandable in biological terms. Simulated data sets can be used for the assessment of novel statistical methods or for the evaluation of the statistical power when designing a study

Distinct Roles of ComK1 and ComK2 in Gene Regulation in Bacillus cereus

The B. subtilis transcriptional factor ComK regulates a set of genes coding for DNA uptake from the environment and for its integration into the genome. In previous work we showed that Bacillus cereus expressing the B. subtilis ComK protein is able to take up DNA and integrate it into its own genome. To extend our knowledge on the effect of B. subtilis ComK overexpression in B. cereus we first determined which genes are significantly altered. Transcriptome analysis showed that only part of the competence gene cluster is significantly upregulated. Two ComK homologues can be identified in B. cereus that differ in their respective homologies to other ComK proteins. ComK1 is most similar, while ComK2 lacks the C-terminal region previously shown to be important for transcription activation by B. subtilis ComK. comK1 and comK2 overexpression and deletion studies using transcriptomics techniques showed that ComK1 enhances and ComK2 decreases expression of the comG operon, when B. subtilis ComK was overexpressed simultaneously

A Bayesian method for evaluating and discovering disease loci associations

Background: A genome-wide association study (GWAS) typically involves examining representative SNPs in individuals from some population. A GWAS data set can concern a million SNPs and may soon concern billions. Researchers investigate the association of each SNP individually with a disease, and it is becoming increasingly commonplace to also analyze multi-SNP associations. Techniques for handling so many hypotheses include the Bonferroni correction and recently developed Bayesian methods. These methods can encounter problems. Most importantly, they are not applicable to a complex multi-locus hypothesis which has several competing hypotheses rather than only a null hypothesis. A method that computes the posterior probability of complex hypotheses is a pressing need. Methodology/Findings: We introduce the Bayesian network posterior probability (BNPP) method which addresses the difficulties. The method represents the relationship between a disease and SNPs using a directed acyclic graph (DAG) model, and computes the likelihood of such models using a Bayesian network scoring criterion. The posterior probability of a hypothesis is computed based on the likelihoods of all competing hypotheses. The BNPP can not only be used to evaluate a hypothesis that has previously been discovered or suspected, but also to discover new disease loci associations. The results of experiments using simulated and real data sets are presented. Our results concerning simulated data sets indicate that the BNPP exhibits both better evaluation and discovery performance than does a p-value based method. For the real data sets, previous findings in the literature are confirmed and additional findings are found. Conclusions/Significance: We conclude that the BNPP resolves a pressing problem by providing a way to compute the posterior probability of complex multi-locus hypotheses. A researcher can use the BNPP to determine the expected utility of investigating a hypothesis further. Furthermore, we conclude that the BNPP is a promising method for discovering disease loci associations. © 2011 Jiang et al

Phenotypic Variation and Bistable Switching in Bacteria

Microbial research generally focuses on clonal populations. However, bacterial cells with identical genotypes frequently display different phenotypes under identical conditions. This microbial cell individuality is receiving increasing attention in the literature because of its impact on cellular differentiation, survival under selective conditions, and the interaction of pathogens with their hosts. It is becoming clear that stochasticity in gene expression in conjunction with the architecture of the gene network that underlies the cellular processes can generate phenotypic variation. An important regulatory mechanism is the so-called positive feedback, in which a system reinforces its own response, for instance by stimulating the production of an activator. Bistability is an interesting and relevant phenomenon, in which two distinct subpopulations of cells showing discrete levels of gene expression coexist in a single culture. In this chapter, we address techniques and approaches used to establish phenotypic variation, and relate three well-characterized examples of bistability to the molecular mechanisms that govern these processes, with a focus on positive feedback.

A Comparative Analysis Shows Morphofunctional Differences between the Rat and Mouse Melanin-Concentrating Hormone Systems

Sub-populations of neurons producing melanin-concentrating hormone (MCH) are characterized by distinct projection patterns, birthdates and CART/NK3 expression in rat. Evidence for such sub-populations has not been reported in other species. However, given that genetically engineered mouse lines are now commonly used as experimental models, a better characterization of the anatomy and morphofunctionnal organization of MCH system in this species is then necessary. Combining multiple immunohistochemistry experiments with in situ hybridization, tract tracing or BrdU injections, evidence supporting the hypothesis that rat and mouse MCH systems are not identical was obtained: sub-populations of MCH neurons also exist in mouse, but their relative abundance is different. Furthermore, divergences in the distribution of MCH axons were observed, in particular in the ventromedial hypothalamus. These differences suggest that rat and mouse MCH neurons are differentially involved in anatomical networks that control feeding and the sleep/wake cycle

Rule based classifier for the analysis of gene-gene and gene-environment interactions in genetic association studies

<p>Abstract</p> <p>Background</p> <p>Several methods have been presented for the analysis of complex interactions between genetic polymorphisms and/or environmental factors. Despite the available methods, there is still a need for alternative methods, because no single method will perform well in all scenarios. The aim of this work was to evaluate the performance of three selected rule based classifier algorithms, RIPPER, RIDOR and PART, for the analysis of genetic association studies.</p> <p>Methods</p> <p>Overall, 42 datasets were simulated with three different case-control models, a varying number of subjects (300, 600), SNPs (500, 1500, 3000) and noise (5%, 10%, 20%). The algorithms were applied to each of the datasets with a set of algorithm-specific settings. Results were further investigated with respect to a) the Model, b) the Rules, and c) the Attribute level. Data analysis was performed using WEKA, SAS and PERL.</p> <p>Results</p> <p>The RIPPER algorithm discovered the true case-control model at least once in >33% of the datasets. The RIDOR and PART algorithm performed poorly for model detection. The RIPPER, RIDOR and PART algorithm discovered the true case-control rules in more than 83%, 83% and 44% of the datasets, respectively. All three algorithms were able to detect the attributes utilized in the respective case-control models in most datasets.</p> <p>Conclusions</p> <p>The current analyses substantiate the utility of rule based classifiers such as RIPPER, RIDOR and PART for the detection of gene-gene/gene-environment interactions in genetic association studies. These classifiers could provide a valuable new method, complementing existing approaches, in the analysis of genetic association studies. The methods provide an advantage in being able to handle both categorical and continuous variable types. Further, because the outputs of the analyses are easy to interpret, the rule based classifier approach could quickly generate testable hypotheses for additional evaluation. Since the algorithms are computationally inexpensive, they may serve as valuable tools for preselection of attributes to be used in more complex, computationally intensive approaches. Whether used in isolation or in conjunction with other tools, rule based classifiers are an important addition to the armamentarium of tools available for analyses of complex genetic association studies.</p

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments.

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the θ_{μe} mixing angle are derived that constitute the most constraining limits to date over five orders of magnitude in the mass-squared splitting Δm_{41}^{2}, excluding the 90%&nbsp;C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90%&nbsp;CL_{s} for Δm_{41}^{2}&lt;13  eV^{2}. Furthermore, the LSND and MiniBooNE 99%&nbsp;C.L. allowed regions are excluded at 99% CL_{s} for Δm_{41}^{2}&lt;1.6 eV^{2}

Examination of polymorphic glutathione S-transferase (GST) genes, tobacco smoking and prostate cancer risk among Men of African Descent: A case-control study

<p>Abstract</p> <p>Background</p> <p>Polymorphisms in <it>glutathione S-transferase </it>(GST) genes may influence response to oxidative stress and modify prostate cancer (PCA) susceptibility. These enzymes generally detoxify endogenous and exogenous agents, but also participate in the activation and inactivation of oxidative metabolites that may contribute to PCA development. Genetic variations within selected <it>GST </it>genes may influence PCA risk following exposure to carcinogen compounds found in cigarette smoke and decreased the ability to detoxify them. Thus, we evaluated the effects of polymorphic <it>GSTs </it>(<it>M1</it>, <it>T1</it>, and <it>P1</it>) alone and combined with cigarette smoking on PCA susceptibility.</p> <p>Methods</p> <p>In order to evaluate the effects of <it>GST </it>polymorphisms in relation to PCA risk, we used TaqMan allelic discrimination assays along with a multi-faceted statistical strategy involving conventional and advanced statistical methodologies (e.g., Multifactor Dimensionality Reduction and Interaction Graphs). Genetic profiles collected from 873 men of African-descent (208 cases and 665 controls) were utilized to systematically evaluate the single and joint modifying effects of <it>GSTM1 </it>and <it>GSTT1 </it>gene deletions, <it>GSTP1 </it>105 Val and cigarette smoking on PCA risk.</p> <p>Results</p> <p>We observed a moderately significant association between risk among men possessing at least one variant <it>GSTP1 </it>105 Val allele (OR = 1.56; 95%CI = 0.95-2.58; p = 0.049), which was confirmed by MDR permutation testing (p = 0.001). We did not observe any significant single gene effects among <it>GSTM1 </it>(OR = 1.08; 95%CI = 0.65-1.82; p = 0.718) and <it>GSTT1 </it>(OR = 1.15; 95%CI = 0.66-2.02; p = 0.622) on PCA risk among all subjects. Although the <it>GSTM1</it>-<it>GSTP1 </it>pairwise combination was selected as the best two factor LR and MDR models (p = 0.01), assessment of the hierarchical entropy graph suggested that the observed synergistic effect was primarily driven by the <it>GSTP1 </it>Val marker. Notably, the <it>GSTM1</it>-<it>GSTP1 </it>axis did not provide additional information gain when compared to either loci alone based on a hierarchical entropy algorithm and graph. Smoking status did not significantly modify the relationship between the <it>GST </it>SNPs and PCA.</p> <p>Conclusion</p> <p>A moderately significant association was observed between PCA risk and men possessing at least one variant <it>GSTP1 </it>105 Val allele (p = 0.049) among men of African descent. We also observed a 2.1-fold increase in PCA risk associated with men possessing the <it>GSTP1 </it>(Val/Val) and <it>GSTM1 </it>(*1/*1 + *1/*0) alleles. MDR analysis validated these findings; detecting <it>GSTP1 </it>105 Val (p = 0.001) as the best single factor for predicting PCA risk. Our findings emphasize the importance of utilizing a combination of traditional and advanced statistical tools to identify and validate single gene and multi-locus interactions in relation to cancer susceptibility.</p

Discovering joint associations between disease and gene pairs with a novel similarity test

Genes in a functional pathway can have complex interactions. A gene might activate or suppress another gene, so it is of interest to test joint associations of gene pairs. To simultaneously detect the joint association between disease and two genes (or two chromosomal regions), we propose a new test with the use of genomic similarities. Our test is designed to detect epistasis in the absence of main effects, main effects in the absence of epistasis, or the presence of both main effects and epistasis. Results: The simulation results show that our similarity test with the matching measure is more powerful than the Pearson's chi(2) test when the disease mutants were introduced at common haplotypes, but is less powerful when the disease mutants were introduced at rare haplotypes. Our similarity tests with the counting measures are more sensitive to marker informativity and linkage disequilibrium patterns, and thus are often inferior to the similarity test with the matching measure and the Pearson 's chi(2) test. Conclusions: In detecting joint associations between disease and gene pairs, our similarity test is a complementary method to the Pearson's chi(2) test

Synergistic Association of PTGS2 and CYP2E1 Genetic Polymorphisms with Lung Cancer Risk in Northeastern Chinese

BACKGROUND: Lung cancer is the most common cause of cancer-related deaths worldwide. The aim of this study was to investigate the association of five extensively-studied polymorphisms in PTGS2 (rs689466, rs5275, rs20417) and CYP2E1 (rs2031920, rs6413432) genes with lung cancer risk in a large northeastern Chinese population. METHODOLOGY/PRINCIPAL FINDINGS: This is a hospital-based case-control study involving 684 patients with lung cancer and 604 cancer-free controls. Genotyping was performed using the PCR-LDR method. Data were analyzed using Haplo.stats and MDR programs. There were significant differences between patients and controls in allele/genotype distributions of rs5275 (P = 0.002/0.003) and rs6413432 (P = 0.037/0.044), as well as in genotype distributions of rs689466 (P = 0.02). The risk for lung cancer associated with the rs5275-C mutant allele was decreased by 60% (95% CI [confidence interval]: 0.21-0.74; P = 0.004) under the recessive model. Carriers of rs689466-G mutant allele had a 28% (95% CI: 0.57-0.92; P = 0.008) reduced risk of developing lung cancer relative to the AA genotype carriers. In haplotype analysis, haplotype G-C-C-T (in order of rs689466, rs5275, rs2031920 and rs6413432) decreased the odds of lung cancer by 28% (95% CI: 0.51-0.93; P = 0.019) after adjusting for confounding factors, whereas haplotype A-T-T-T had 1.49-fold (95% CI: 1.21-1.79; P = 0.012) increased risk for lung cancer. Using MDR method, the overall best model including rs5275, rs689466 and rs6413432 polymorphisms was identified with a maximal testing accuracy of 66.1% and a maximal cross-validation consistency of 10 out of 10 (P = 0.003). CONCLUSIONS/SIGNIFICANCE: Our findings demonstrated a potentially synergistic association of PTGS2 and CYP2E1 polymorphisms with the underlying cause of lung cancer in northeastern Chinese